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Abstract. Exact semianalytic solutions for GRB outflows were recently derived using the equations 



\ of special-relativistic ideal MHD (see the contribution by Vlahakis & Konigl for a summary). This 



contribution focuses on the implications of these results to various modeling and observational 
issues in GRB sources, including the baryon loading problem, polarization measurements of the 
prompt and reverse-shock emission, and the possible existence of a two-component outflow. 
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^ : 1. MAGNETIC DRIVING OF GRB OUTFLOWS 
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Gamma-ray burst (GRB) outflows are likely powered by the extraction of rotational en- 
ergy from a newly formed stellar-mass black hole or neutron star, or from a surrounding 
debris disk established in the course of the central object's formation [e.g.,[l5|]. Magnetic 
fields threading the central object or disk provide the most plausible means of extract- 
ing the inferred amount of energy on the timescale of the burst; they can also guide, 
collimate, and accelerate the flow [see 20, and references therein]. This picture is sup- 
ported by a recent measurement of a high (80 ± 20%) linear polarization in the prompt 
y-ray emission from GRB 021206 [5], which can be plausibly interpreted in terms of a 
large-scale magnetic field advected from the origin [e.g., lull and is consistent with mag- 
netic driving by an ordered field that threads the source. Although purely hydrodynamic 
driving powered by neutrino emission or magnetic field dissipation at the source can 
probably be ruled out [e.g.,[a,lu], thermal effects may nonetheless dominate the initial 
acceleration of magnetic jets re.g.. ll6l 12011 . 1 

Motivated by the above considerations, Vlahakis & Konigl ll2ll I22L hereafter VK03a 
and VK03b, respectively] constructed a general formalism for special-relativistic ideal 
MHD, allowing for the presence of a baryonic component as well as of a "hot" electron- 
positron/radiation component that can dominate the pressure. They showed how one can 
derive exact semianalytic solutions for axisymmetric outflows under the assumption of 
radial self-similarity and presented illustrative results for representative GRB parame- 
ters. Vlahakis, Peng, & Konigl H23L hereafter VPK03] further generalized this scheme 
by obtaining solutions for initially neutron-rich outflows, which they used to address the 
baryon loading problem in GRB source models. A general description of the formal- 



1 It has also been argued [e.g., 9] that electromagnetic energy dissipation could contribute to the conver- 
sion of Poynting energy into kinetic energy throughout the acceleration region of such flows. 



ism and of the derived solutions is given in Vlahakis & Konigl's contribution in these 
Proceedings. The present contribution provides a brief overview of the main results and 
focuses on their observational implications. 



2. RELATIVISTIC MHD SOLUTIONS 

The initial (subscript i) magnetic field amplitude can be inferred from an estimate of the 
injected energy, S\ = (Poynting flux) x (surface area) x (burst duration). In a disk geom- 
etry [with initial cylindrical radius GJj and radial width (ACJ)j], S\ ~ cE\B§ jGJ^AQJ^A?, 
where the electric field is given by E = B V V^/ c — B^Vp/ c (with the subscripts p and 
denoting the poloidal and azimuthal components, respectively). For characteristic 
parameter values \S\ ~ 10 52 ergs, G5j ~ (AC7)i ~ 10 6 cm, At ~ 10 s], one obtains 
5j ~ 10 14 — 10 15 G. This field is most plausibly generated by differential-rotation am- 
plification of a much weaker poloidal field component that originally threads the source. 

If |5 p j/5^ ;i | > 1, a trans-Alfvenic outflow is produced, whereas if \B$^/B V j\ > 1, 
the outflow is super- Alfvenic from the start. The latter situation may correspond to 
amplified toroidal flux loops that have been disconnected by magnetic reconnection 
and escape from the disk surface in a nonsteady fashion. Exact solutions for these two 
situations were derived in VK03a and VK03b, respectively. It was demonstrated that, 
in either case, Poynting flux-dominated jets can transform ^ 50% of their magnetic 
energy into baryon kinetic energy (with £k ~ 10 51 ergs and terminal Lorentz factors 
7oo ~ 10 2 — 10 3 ). If relativistic pairs and radiation dominate the initial enthalpy, 
then a thermal acceleration zone develops at the base of the flow and remains dominant 
until the specific enthalpy drops below ~ c 2 , at which point magnetic acceleration takes 
over. In contrast to the trans-Alfvenic solutions, part of the enthalpy flux in the super- 
Alfvenic flows is transformed into Poynting flux during the thermal acceleration phase. 
Furthermore, the subsequent, magnetically dominated acceleration in these flows can be 
significantly less rapid than in the trans-Alfvenic case. 

The derived solutions have a free parameter, F, which controls the distribution of 
the poloidal current / = cWB^/l. For F > 1 the flow is in the current-carrying regime, 
with the poloidal current density being antiparallel to the magnetic field. In this case 
the current tends to zero as the symmetry axis is approached, so such solutions should 
provide a good representation of the conditions near the axis of a highly collimated flow. 
Conversely, solutions with F < 1 correspond to the return-current regime (in which the 
poloidal current density is parallel to the field) and are most suitable at larger cylindrical 
distances. Although the detailed global current distribution cannot be modeled using the 
self- similarity approach, one can nevertheless generate "hybrid" flow configurations that 
combine a current-carrying solution for low values of U5 and a return-current solution for 
high values of UJ (see Fig. 1 below for an example). Initially Poynting-dominated flows 
that attain a rough equipartition between the kinetic and Poynting energy fluxes at large 
distances from the origin have F close to 1. When F > 1 the Lorentz force can collimate 
the flow to cylindrical asymptotics. For F < 1 the collimation is weaker and the flow 
only reaches conical asymptotics; however, the acceleration is more efficient in this case 
in that a larger fraction of the Poynting flux is converted into kinetic energy. 



3. IMPLICATIONS TO THE BARYON LOADING PROBLEM 



As an illustration of the unique properties of the relativistic MHD solutions, consider 
the ramifications of a hydromagnetic jet model to the baryon loading problem in GRB 
outflows. The apparent difficulty stems from a comparison between the estimated mass 
of protons in the jet, M pro ton = 3 x 10~ 6 (E K /10 51 ergs)(7oc/200) _1 M , and the mini- 
mum mass of the debris disk from which the jet is thought to originate, obtained under 
the assumption that at most ~ 10% of the disk gravitational potential energy could be 
converted into outflow kinetic energy. This comparison implies that the outflow can 
comprise at most ~ 10~ 4 of the disk mass, whereas disk outflow models that utilize a 
large fraction of the disk potential energy typically also entail substantial mass loading. 
One approach to this issue has been to postulate that the outflow emerges along mag- 
netic field lines that thread the black-hole event horizon and not the disk, but then the 
converse problem — how to avoid having too few baryons — must be addressed [e.g., 
fl3ll . A possible resolution of the problem in the context of disk- fed jet models was pro- 
posed in [10], where it was noted that such outflows are expected to be neutron-rich 
[neutron/proton ratios as high as n/p ~ 20 — 30; e.g.,[Tsll3,|2j]. Since only the charged 
outflow component couples to the electromagnetic field, the neutrons could potentially 
decouple from the protons before the latter attain their terminal Lorentz factor. In this 
picture, the inferred value of M proton may represent only a small fraction of the total 
baryonic mass ejected from the disk, which would alleviate the loading problem. How- 
ever, it can be shown that, for purely hydrodynamic outflows, the Lorentz factor y& at 
decoupling is at least a few times 10 2 [e.g.,|7|,|2j,|23|]. This implies that /d/y^ ~ 1 and 
hence that the protons end up with only a small fraction of the injected energy, which is 
not a satisfactory resolution of the problem. 

As demonstrated by VPK03, the incorporation of magnetic fields makes it possible 
to attain yj <C y^ and thereby reclaim the promise of the Fuller et al. proposal. They 
wrote down the equations of motion for the neutron component (which couples to 
the protons through a collisional drag) and for the charged component (incorporating 
protons and their neutralizing electrons as well as initially "hot" pairs and radiation), 
and simplified them by considering a well-coupled neutral/charged fluid for y < 7d and 
only the charged fluid component for y > y^. The pre-decoupling region was described 
by a super- Alfvenic outflow solution. As noted in § 2, in this case part of the enthalpy 
flux is converted into Poynting flux during the initial thermal acceleration phase. This 
reduces the acceleration rate, so at the point of decoupling (when Vp ro ton — ^neutron ~ c) 
the Lorentz factor is still comparatively low. The energy deposited into the Poynting 
flux is returned to the matter beyond the decoupling point as kinetic energy, thereby 
enhancing the acceleration efficiency of the proton component. The end result is a 
large y^/ yj ratio and comparable terminal kinetic energies in the proton and neutron 
components, in clear contradistinction to the purely hydrodynamic solutions. 

An illustrative solution is shown in Fig. I. 2 The top panel shows the behavior of the 



In this example n/p = 30, the pre-decoupling and post-decoupling regions correspond to the current- 
carrying (F = 1.05) and return-current {F = 0.1) regimes, respectively, and the flow collimates from an 
initial opening half-angle of 55° to 0j « 20°. 
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FIGURE 1. Illustrative relativistic-MHD solution of a neutron-rich outflow, (a) Components of the total 
energy flux, normalized by the mass flux x c 2 , as functions of height along a fiducial magnetic field line. 
The Poynting and enthalpy curves are discontinuous at the decoupling point, reflecting the decrease in the 
mass flux of field-coupled gas above that point, (b) Components of the proton-neutron drift velocity. 



various components of the energy flux, corroborating the qualitative description given 
above. The thermal acceleration effectively terminates at a height z ~ 10 9 cm above 
the disk, and the neutrons decouple from the protons at Zd ~ 10 13 cm, corresponding to 
7d ~ 15. By the time of decoupling the neutrons have acquired ~ 2/3 of the injected 
energy, with the remainder residing predominantly in the electromagnetic field. The 
latter portion is then transferred with almost 100% efficiency into proton kinetic energy, 
so that, ultimately, the protons have = 200 and £k,proton ~ 10 51 ergs 0.5 /^neutron- 
The proton jet thus carries ~ 1/3 of the injected energy but only ~ 3% of the injected 
mass. The lower panel of Fig. 1 shows that, even though the decoupling in this case is 
initiated by the growth of the n-p drift velocity along the poloidal magnetic field, there is 
also a transverse drift component (induced by the ongoing magnetic collimation), which 
at the time of decoupling is V n eutron, _i_ ~ 0. 1 c. 3 



3 The exact value of the angle between V n and V p at decoupling can only be obtained by solving the 
equations of motion without the "strong coupling" approximation adopted in the solution shown in Fig. 1. 





4. ADDITIONAL IMPLICATIONS 



Polarization — Magnetic driving of GRB outflows by large-scale, ordered magnetic 
fields would naturally lead to a large linear polarization P in the prompt y-ray emission, 
and a high value of P is also predicted for the emission from the reverse shock (the 
"optical flash" and "radio flare") 111 211 . As shown in 111 ill , in this picture the prompt emis- 
sion may be expected to exhibit P ~ 43% — 61% for typical values of the synchrotron- 
radiation spectral index, consistent with the observations of GRB 021206 a]. The or- 
dered field is also expected to induce measurable circular polarization ifbUl . 
Two-Component Outflow — The decoupled neutrons in a neutron-rich outflow will 
undergo /3 decay into protons at a distance ~ 4 x lO 1 ^ (W 15) cm. In contrast with 
the situation in purely hydrodynamic outflow models iTlvL Ql], there may well be no 
interaction between the two decoupled components in the MHD case since their motions 
are not collinear (see Fig. lb). The latter scenario thus gives rise to a 2-component 
outflow: an outer (wider) component (comprising the decoupled neutrons) that carries 
most of the energy and may be responsible (after the neutrons decay) for the bulk of 
the optical/radio afterglow, and an inner (narrower) component (comprising the original 
protons) that accounts for the prompt y-rays and possibly also for much of the X-ray 
afterglow. A 2-component outflow of this type was inferred in GRB 030329 (3lll9J]. A 
more detailed investigation of this scenario is currently under way. If £k,narrow ^ £k,wide 
and 0j, narrow /^j, wide ^ 1/3, this picture would make it possible to reconcile current 
inferences of the radiated y-ray energy [e.g., 4] with internal-shock models. 
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